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When  new  design  criteria  are  developed  it  is  often  necessary  to 

develop  new  information  which  will  ultimately  be  tabulated  for  the  practicing 

engineer  who  wishes  to  avoid  the  tedious  calculations  required  to  select 

optimum  sections.  In  this  paper  we  will  investigate  the  possibility  of 

using  geometric  programming  methods  to  optimize  the  geometry  of  a pitched 

tapered  beam  subjected  to  uniformly  distributed  loads.  The  beam  geometry 

is  shown  in  Figure  1 and  is  characterized  by  the  span,  2L,  the  width  b, 

the  radius  of  curvature  R,  the  heights  at  the  center  H and  Hc,  the 

height  at  the  support  H , the  roof  angle  P,  and  the  slope  of  the  lower 

s 

surface  cp. 

In  this  paper  we  have  chosen  to  optimize  the  design  of  the  pitched 
tapered  beam  with  respect  to  volume.  Fox  [1]  presented  a computer  program 
for  minimizing  the  volume  of  beams  using  a technique  that  required  a 
detailed  understanding  of  the  constraints  and  employed  a rather  arbitrary 
fixed-step  search  technique  to  move  along  constraint  surfaces  until  a 
minimum  volume  was  achieved.  This  method  provided  optimal  designs  in  many 
cases  but  because  of  the  arbitrary  step  sizes  used  in  the  program  it  can 
be  shown  that  the  attaining  of  minimum  volumes  is  not  always  ensured. 

Since  this  program  was  developed  it  has  been  shown  by  Barrett  et  al.  [2] 
that  the  tensile  strength  perpendicular  to  grain  of  timber  is  not  inde- 
pendent of  member  geometry  as  has  been  previously  assumed.  In  particular, 
given  a set  of  design  parameters  x^  which  characterize  the  beam  geometry, 
the  stresses  and  deflection  are  checked  using  formulae  normally  specified 
in  building  codes.  Stresses  and  deflections  so  calculated  must  not  exceed 
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the  maximum  (or  allowable)  values.  These  relationships  can,  therefore, 
be  expressed  in  the  form  of  inequality  constraints 


(1.1) 


where  cr^  is  the  computed  stress  component  and  is  the  maximum 

allowable  value  for  that  stress  component.  In  addition  to  stress  con- 
straints, deflection  constraints  are  usually  employed  to  keep  the 
deflection  D of  the  member  below  a specified  fraction  of  the  span  and 
accordingly  the  constraint  will  have  the  form 


D(xjL)  < aS 


(1.2) 


where  a is  the  appropriate  fraction  and  S is  the  beam  span. 

In  general,  there  will  be  other  relations  which  control  beam 
geometry.  For  example  in  beam  design  often  the  ratio  of  beam  height  to 
width  is  restricted  to  satisfy  lateral  stability  requirements.  In  the 
case  of  the  pitched  tapered  beam  there  will  be  requirements  for  a minimum 
length  of  tangent  on  the  beam  lower  surface.  These  constraints  can  be 
written  in  the  form 


W £ 1 • 


(1.5) 


The  allowable  tensile  strength  perpendicular  to  grain  for  pitched 
tapered  beams  is  given  by 
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F = 

r 


(1.4) 


where  h is  the  beam  width,  H the  height  at  the  centerline,  R the 

m 

radius  to  midheight  and  cp  the  angle  in  radians  between  the  centerline 
and  the  tangent  point.  The  allowable  stress  according  to  Eq.  (1.4)  is 
generally  lower  than  the  65  psi  used  previously,  thereby  requiring 
modification  of  the  beam  configurations  which  have  been  tabulated  in  [3). 

The  optimal  beam  geometry  will  be  found  by  formulating  the 
optimization  problem  as  a generalized  geometric  (signomial)  programming 
problem.  In  the  next  section  the  detailed  engineering  formulation  of 
the  problem  will  be  presented.  The  signomial  programming  formulation  is 
derived  in  Section  3 and  in  Section  4 the  results  of  a few  sample 
designs  are  given. 
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2.  Engineering  Considerations 

Given  specifications  for  the  uniformly  distributed  design  dead 
load  and  snow  load  oig,  the  total  span  2L  and  the  roof  angle  3, 

the  design  engineer  must  develop  a beam  of  adequate  capacity.  The 

I1  capacity  of  a beam  is  assessed  by  evaluating  stress,  deflection  and  geometric 

h 

criteria  which  are  normally  specified  by  building  code  authorities  or 
dictated  by  manufacturing  requirements.  For  this  paper  the  design 
criteria  for  allowable  shear  and  bending  stress,  F and  F_  and  the 
deflection  constraints  are  those  specified  by  the  Canadian  Standard 
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Association  Standard  CSA  086-1970  [U],  Currently  the  allowable  tension 
perpendicular  to  grain  stress  F , permitted  in  Canada  for  Douglas  fir, 
dry  service  condition  and  normal  load  duration  is  65  psi.  It  has  been 
shown  that  this  allowable  stress  is  not  adequately  conservative  for  large 
beams  and  for  this  paper  the  allowable  stress  shall  be  that  recommended 
by  Barrett  et  al.  [2]. 

In  design  the  engineer  strives  to  produce  a minimum  cost  structure 
consistent  with  the  requirements  for  public  safety  and  protection  of 
property.  This  may  be  accomplished  by  minimizing  the  volume  of  material 
used  in  the  members.  Accordingly,  in  our  optimization  problem  the 
objective  function  to  be  minimized  is  the  volume  of  the  beam  subject  to 
the  design  constraints  developed  below. 

There  are  constraints  on  the  bending  stresses  at  the  beam  center- 
line  and  tangent  point  (TP,  Figure  1),  the  shear  stress  at  the  support 
(A,  Figure  1)  and  tension  perpendicular  to  grain  stress  at  the  centerline. 
The  specific  form  of  these  constraints  for  the  pitched  tapered  beam 
subjected  to  a uniformly  distributed  load  is  as  follows. 

For  the  bending  stresses  ct-q,  we  require  cr_  < F_  where  F_  is 
the  given  allowable  bending  stress.  Bending  stresses  must  be  checked  at 
two  positions,  at  the  centerline  and  at  the  tangent  point.  Therefore  at 
the  centerline  the  constraint  has  the  following  form 

6m 

Cg  = (1  + 2,7  tang]  < Fg  (2.1) 
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where  Mc,  the  bending  moment  at  the  centerline  is  given  by 


ali 


Mc  = 2T 


(2.2) 


Similarly,  the  bending  stress  at  the  tangent  point  is  given  by 


6Hn 

— i [1  + 2«7  tan  p]  < F 

*4  E 


(2.3) 


where  Hp,  the  bending  moment  at  the  tangent  point  is  given  by 


Mj  = Mc  - <oR2  sin2  q>/2l+  . 


(2.1+) 


Here, 

H,  H^,  = beam  heights  at  centerline  and  tangent  point,  respectively, 
inches, 

L = half-span,  inches, 
b = beam  width,  inches, 
p = roof  angle,  degrees, 

od  = uniformly  distributed  load,  pounds  per  lineal  foot. 


The  beam  loading  tends  to  increase  the  radius  of  curvature  of  the 
beam,  thereby  introducing  stresses  in  the  radial  direction.  The  magnitude 
of  the  tension  perpendicular  to  grain  stress  cr  is  given  by 


6 


where 


and 


a 

r 


6m 


K = A + B(#-) 


m 


( H / 

c(r) 

m 


(2.5) 


(2.6) 


R = radius  to  midheight  (R  + H/2) , inches 
m 

A,B,C  = functions  of  3,  tabulated  in  [4]. 

The  shear  stresses  t at  the  support  are  required  to  satisfy  the 

relation  t < F and  this  constraint  is  formulated  as  follows 
— s 


where 


and 


T = l -X_  < F 
2 bH  - s 
s 


(2.7) 


(2.8) 


H = beam  height  at  the  support,  inches 
s 

v = shear  force,  pounds 


A constraint  on  the  midspan  deflection  of  beams  is  normally  imposed 

to  prevent  excessive  deflection  which  could  damage  ceiling  materials.  The 

maximum  allowable  deflection  6 , is  usually  expressed  as  a function  of 

max 

the  total  span.  The  corresponding  constraint  in  our  case  is 

k 

* _ _ 

288EI  - max 

where 

Y = 0.2  + 0.6  H /H  , (see  [4])  (2.10) 

c s 

I = bHI * 3/l2  (2.11) 

c 

and 

E = modulus  of  elasticity,  psi 

Hc  = height  at  centerline  for  the  double -tape red  component  of 


the  beam,  inches. 


*> 

l 

t The  final  constraints  are  constraints  on  geometry.  The  first  is 

| introduced  to  ensure  that  the  tangent  point  (T.P.,  Figure  l)  is 

i 

< positioned  at  an  adequate  distance  from  the  end  of  the  beam.  This 

constraint  is  required  to  prevent  springback  when  the  beam  is  released 
from  the  clamps  after  curing  of  the  glue  is  complete.  For  convenience 
the  constraint  is  formulated  as  follows 


L - R sin  cp  > qH 


(2.12) 


where  a is  a given  constant. 

To  complete  the  engineering  formulation  it  is  necessary  to  specify 
the  allowable  values  to  be  used  in  the  stress  and  deflection  constraints. 
These  values  depend  on  the  species  of  wood  used  in  the  beam.  If 
Douglas  fir  is  used  the  allowable  values  are  as  follows 


Fd  = 2760(1  - 2000 (t/R)2)  (psi ) 

D 


where 


t = lamination  thickness,  inches. 


F = K(bHF.  Cp) 

r m 


-0.2 


F = VjO 
s 


where 


K = 324  (if  ■'(  is  in  radians) 


8 = 

max  5 


where  we  assume  ? = 180. 


(psi ) 
(psi  ) 


(2.15) 


(2.14) 

(2.15) 


[2.16) 
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The  above  formulas  for  F and  F include  a 15 $ increase  in  allowable 

B s 

stress  for  snow  load  conditions. 

Finally,  the  radius  of  curvature  R is  constrained  to  be  greater 
than  or  equal  to  330  inches  so  that  excessive  stresses  will  not  be 
introduced  when  the  beams  are  fabricated. 

The  variable  portion  of  the  half -beam  volume,  to  be  minimized,  is 
given  by  (see  Figure  l) 

V = Lb  (H  + H ) + bR^ (tan  <p  - <p)  . (2.17) 

s c 

2 

The  total  half -beam  volume  is  equal  to  V - L b tan 

This  concludes  the  engineering  formulation  of  the  optimal  design 
problem. 


3-  Signomial  Programming  Formulation 

In  this  section  we  present  a signomial  programming  formulation  of 
the  optimal  beam  design  problem.  The  main  difficulty  in  reformulating 
this  problem  as  a signomial  program  is  that  some  of  the  relations  appear- 
ing in  the  preceding  section  are  equations,  whereas  signomial  programming 
constraints  must  be  inequalities.  As  will  be  shown  below,  some  of  the 
equations  of  the  engineering  formulation  are  used  to  eliminate  variables 
and  others  are  converted  into  inequalities  in  such  a way  that  hopeful?  • 
they  hold  as  equations  at  an  optimum.  Another  minor  problem  arises  from 
the  appearance  of  trigonometric  functions  in  the  design  formulas.  Taylor 
series  approximations  of  these  functions  are  used  below  to  obtain 
generalized  polynomials  as  required. 

y 


Note  that  converting  an  equality  constraint  g(x)  - 1 into  two 
inequalities  g(x)  < 1,  g(x)  >lis  not  practical,  since  the  algorithm 
used  for  the  numerical  solution  is  based  on  the  assumption  that  the 
interior  of  the  constraint  set  is  nonempty.  For  this  reason,  equalities 
can  only  be  converted  into  one-sided  inequalities.  The  sense  of  the 
inequalities  is  usually  determined  by  physical  or  design  considerations 
(see,  for  example,  [5])*  Unfortunately,  there  are  cases  where  these 
considerations  are  quite  complex  and  cannot  be  observed  by  simple 
inspection  of  the  constraints.  Consequently,  a trial-and-error  approach 
is  necessary.  Simple  examples  can,  however,  be  constructed  showing  that 
not  every  equality  constrained  problem  can  be  solved  by  converting 


equations  into  single  inequalities. 

Let  us  derive  now  the  signomial  program  for  the  beam  design  in 


detail.  The  volume  of  the  beam  as  given  in  the  engineering  formulation 


V = Lb  ( H +•  H ) + blf  (tan  % - Cp)  . 
sc 


O.l) 


Using  the  identity 


H *-  L tan  cp  = H + L tan  8 
c s ^ 


0-2) 


and  a three-term  Taylor  expansion  of  tan  cp 


13  2 5 

tan  cp  = cp  + - cjr  + — cp 


we  obtain  the  volume  to  be  minimized 


V = 2LbH  - L2b  tan  8 +■  L2bcp  + i L2bq r 
c 3 


225  1 23  2 25 

+ ~ LOV  + \ bR  eg  r f- 


10 


where  H , the  beam  height  at  the  support,  has  been  eliminated. 
The  variable  H is  defined  by  the  identity 


H - R 


cos  cp 


0.5) 


In  converting  this  relation  into  an  inequality  we  can  ensure  that  in  an 
optimal  solution  the  inequality  will  hold  as  an  equation  by  writing 


RU  - 


0.6) 


Since  in  (3*4)  we  try  to  lower  the  value  cf  li  as  much  as  possible, 
the  inequality  in  (J.6)  will  be  tight  in  an  optimal  solution.  Substitute 


cos  cp 


- 1 - ~ <P" 


5 4 . f.  1 4 

rr-  Cp  ;rr—  Cf 
24  I 20  H 


0.7) 


intc  O.'  1 and  rearranging  yields 


H •-  Rr/  i-  ~~  RfA  *-  -ii  t.'0 

c 2 4 ?>+  * 720  4 


, 2,  U) 

V .?  • ^ ’ 


,-l  -2.,  r n-l  -2„  10  2 102  4 

i 0 n-  -rP  q-  IIc  - .pr-  9 - ~ — q < 1 


The  bending  stress  onstrain'  a*  centerline  is  given  by  If . l ) , 
and  (2 . 15  ) as 


6m  vl  P.Ttan  0) 

— — — < 276o[  1 - 2000  {\f] 


.3.10  ) 


O O r-\  s~) 

FT.  H r_  +•  OOCOOr  ' • 1 


where 
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(1  + 2.7  tan  e)o- 


4 V27G0  )b 


(5.12) 


Turning  now  to  the  bending  stress  constraint  at  tangent  point  we  have 
from  (2.2),  (2.3),  (2.4)  and  (2.13) 


, (JDL2  <AR2  sin2tp\  (1  + 2.7  tan  p) 
24  ' 24  b£ 


2760[  1 - 2000  (|)"]  >3.13  ) 


2-2  ? 2 -2 
PL  ” - TR  sin  rp  H^,'  + 


2000t2R'c  1 . 


0-14! 


The  beam  height  at  tangent  point  H , is  related  to  the  other  beam  heights 


H,  - (1  - - s-7n-JK)  H + (5_®iS_S:)  H . 
T L c L s 


(3.15  » 


From  (3.2)  and  (3.3)  we  obtain 


ii, , - H + R - R cos  cp  - R tan  p sin  0 . 


<,3.  It ) 


Instead  of  converting  (3.16)  into  an  inequality  and  guessing  its  sense, 
we  substitute  (3-16)  into  (3.14)  and  by  letting 


13  15 

sin  rp  - cp  - qr  + Y2Q  V 


i3. 17) 


12  14 

-os  Ip  =-  1 - - cp  + ^ q> 


(3.18) 


we  write  the  bending  stress  constraint  at  tangent  point  as 


12 


TL  K + H Kcp  + 2 tan  p H Rep  + —f i-f-  H Ttep'  + 2000t  E 


+ 500t2  tan  p H_2cp5  + 2v--— tan  K_2R2cpD  + tan  p H “R' 


-®"V 


+ 2000t2  H~1R"1cp2  + 200Gt ^ -tan  P H-1R"\p^  +■  tan2  p 2000t2H_2Cp2 


(3A  - tan2  p)  2000t  „-2  U _ H-l_  ? tan  6 U-1  * 


H (?)  - H Rep  - -r-^  H Rep' 

- 2C00t2  tan  pH^cp5  - (r  + tan2p)  H~^R2qp2  - (-—-"  T ^ tan — &)  H~2R2rr 


tan  p . -2r2  5 2000t2  .,-1-1  k , nrsnJ2  . „ ,-L-l 

- — p— ^ H R ep  - H R <p  - 4000t  tan  p H 1 ' 


CP 


_ 2000V_tan__£  . 8J1  - tan^.p'  g0OCtg.  ^ 6 ^ ^ 


IT 


Next  we  consider  the  constraint  on  tension  perpendicular  to 
grain  stress.  We  have  by  (,2.5),  12. t)  and  (2.1U) 


6M  A -M  B 6M  C 

h 2 +■  H iR  1 + -A-  R 

b b m b m 


U>HR  c?) 
m 


0.2 


(5.20) 


or 


6M=A~  v- v-2  - 24  h-o-8b-°-v-2  - 24-  < i 

m ^ 0.cT,  m ,0.o„  m - 

D K 


b°‘wK 


where 


R = R + =■  H 
m 2 


(3.21) 

(3-22) 


Direct  substitution  of  (3>2P)  into  (3.21 ) would  yield  a nonsignomial 
constraint  and,  therefore,  should  be  avoided.  Consequently,  we  must  treat 
(3.22)  as  an  inequality  that  has  to  hold  as  an  equation  at  optimum. 


-^W 
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l 

i- 

'(■ 


I 

f| 

rf 


! 


Since  has  both  positive  and  negative  exponents  in  (3-21)  uhe 

sense  of  the  inequality  to  replace  (3.22)  cannot  be  determined  in  advance. 
Multiplying,  however,  the  left  hand  side  of  (3.2l)  by  the  identity 
(R  +•  h/2)/e^  yields  the  new  tension-perpendicular-to  grain  stress  con- 
straint 


6m  A 
c 


,-1.8_-0.c  0.2 


6M  b 


o73Z  RH  Rm  * • ,0.8. 


b K 


b K 


„„-G.8  -1.8  0.2 
RH  R CD 
m 


. 6V  eh0.2b-2.8  0.2  . „-0.8k-0.8  0.2 

0.8„  m .0.8,,  m 


b K 


b K 


3McB  0.2  -1.8  0.2  ?McC  1.2-2.80.2  . 

+■  7TT-  H R <p  t , 7?-  H R__  C-  < 1 

b°.SK  m b0.8K 


m 


(3.23) 


where  all  the  terms  on  the  left  hand  side  are  positive  and  the  exponents 
of  Rm  are  all  negative.  Now  we  convert  (3.22)  into  the  inequality 

Rm  - R + 2 H (3-24) 

or 

R R-1  - r HR-1  < 1 . v*>.25 ) 

m2  — 

Note  that  the  above  considerations  are  valid  if  the  inequality  (3.23)  is 
tight  in  the  optimal  solution.  It  may  happen,  however,  that  both  (3.23) 
and  (3.25)  are  strict  inequalities  at  optimum.  In  this  case  the  sense 
of  (3.25)  must  be  reversed  (such  a reversal  was  in  fact  necessary  in  one 
of  the  cases  solved). 

The  shear  stress  constraint  is  formulated  from  (2.7)  by  using 
(3.5)  and  (3.3).  We  obtain 


a 


I 


a 


i 


? 
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+ L tan  p)H  - I4H  <j>  H < 1 

8bF  hc  ^c  3 c 15  c - 


(3.26: 


The  deflection  constraint  is  formulated  from  12.9),  (2.10),  (2.1l)  and 


(2.16)  as 


T)H~^  + l4.T)lf2H-1  < 1 
* c c s — 


(3.27) 


0.208gmL^ 


(3.28; 


Multiplying  both  sides  of  (3-27)  by  H and  substituting  (3*2)  yields 

s 

iH  2 - t|L  tan  3 H ^ + iL  tan  9 H ^ + Ini  H 2 


< H - L tan  B + L tan  cp 
— c 


(3*29) 


and  by  (3.3)  we  obtain 


51  H~3  + iLcpIi'11"  + j llVtT^  + r,  Lq?lfk  + L tan  3 H~'L 

- 1 L tan  3 - LcpH"1  - j LcpV1  - Lcp5^1  < 1 . (3-30) 


The  constraint  on  the  geometry  of  the  beam  given  by  (2.12)  is  rearranged  to 


? H + 7 B sin  cp  < 1 

L L — 


(3.31) 


and  by  (3-17)  it  becomes 


The  last  geometry  constraint,  R > 330, is  not  treated  explicitly,  since 
the  numerical  algorithm  for  the  solution  of  the  beam  design  problem 
requires  upper  and  lower  bounds  on  all  the  design  variables,  thus  the 
value  330  will  be  used  as  the  lower  bound  on  R. 

The  optimal  beam  design  is  obtained  (after  specifying  the 
appropriate  constants)  by  solving  the  signomial  program  of  minimizing 
(3-4)  subject  to  the  constraints  (3-9),  (3-ll),  (,3.1b),  (3.23),  ('.25), 
(3-26),  (3.30)  and  (3-32).  The  variables  to  be  determined  by  the 
optimization  are  H,  H , R,  R , V and  <p.  Note  that  an  optimal  solution 
to  the  signomial  programming  formulation  of  the  design  problem  is 
acceptable  only  if  the  inequalities  (3. 9), and  (3-25)  hold  as  equations 
at  optimum. 

A few  sample  problems  of  optimal  beam  designs  were  solved  by 
the  computer  code  GGP,  based  on  the  generalized  geometric  (signomial) 
programming  algorithm  of  Avriel,  Dembo  and  Passy  [5].  These  optimal 
design  solutions  are  presented  in  the  next  section. 


4.  Sample  Designs 


Optimal  beam  configurations  are  sought  for  three  different 
spans  and  loading  conditions.  The  specified  beam  parameters  are  shown 


in  Table  1. 


Table  1 


Case 

Input  Parameters 

for  Beam 

Optimization 

Load 

to  (lb/ft  i 

Roof  Angle 
3 (degrees ) 

Half-Span 
L (inches) 

Lamination 
Width  Thickness 

b (inches)  t (inches) 

1 

9.46 

360 

6. 75 

1.5 

1200 

2 

9.46 

240 

6.75 

1.5 

1200 

3 

9.46 

120 

3.00 

1.5 

400 

For  the  above  roof  angle  the  corresponding  constants  are  A --  0.0367, 

B = 0. 

0794,  C = 0.213.  In  addition,  the 

modulus  cf  elasticity  is  assumed 

to  be 

E = 1.93  x 

10 ^ psi  and 

a value  of 

a = 1.5  is  taken  in  (2.12). 

Optimal  solutions  were  obtained  by  the  computer  code  GGP  in  less 
than  10  seconds  of  CPU  time  on  an  IBM  37o/l66  computer.  The  optimal 
design  variables  are  listed  in  Table  2. 


Table  2 

Optimal  Design  Variables 


Case 

Volume 

ft3 

9 

degrees 

H 

inches 

«c 

inches 

R 

inches 

inches 

1 

184.32 

3.64 

70.7 

68.3 

1063 

1099 

2 

71.05 

5-01 

50.3 

47.0 

860 

866 

3 

5.93 

6. 16 

19-5 

17.6 

330 

340 
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It  is  interesting  to  observe  the  binding  design  constraints  at  optimum 
for  the  above  cases  (in  addition  to  those  which  must  be  tight  because 
they  were  originally  equations). 


Table  3 

Binding  Design  Constraints  at  Optimum 

Case  Binding  Design  Constraints 

1 Tension  -4  grain  stress  (3*23);  Shear  stress  (3.26) 

2 Tension  -L  grain  stress  (3-23);  Shear  stress  (3 .26) 

3 Bending  at  tangent  point  (3-19);  Shear  stress  (3.26) 


In  Cases  1 and  2 constraint  (3-23)  is  tight  at  optimum  and  consequently 
(3.25  ),  the  defining  relation  for  R^,  is  also  satisfied  as  an  equation. 

In  Case  3>  however,  (3-23)  is  no  longer  binding  and  at  first  we  obtained 
a solution  in  which  both  (3*23)  and  (3>25)  were  strict  inequalities.  We, 
therefore,  reversed  the  sense  of  the  inequality  in  (3-24)  and  (3-25)  to 

R > R + | H (4.1) 

m — 2 

and 

RR"1  + i HR"1  < 1 (4.2) 

m 2 m- 

respectively,  and  (3-25)  was  replaced  by  (4.2)  in  the  program.  This  change 
resulted  in  the  above  listed  optimal  solution  for  Case  3 in  which,  of 
course,  (4.2)  held  as  an  equation. 


j 

1 


1 
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